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^ ^ ' Abstract. Globally regular (ie. asymptotically flat and regular interior), spher- 

^ ^ ically symmetric and localised ( "particle-like" ) solutions of the coupled Ein- 

■ stein Yang-Mills (EYM) equations with gauge group SU{2) have been known 

for more than 20 years, yet their properties are still not well understood. Spher- 
ically symmetric Yang-Mills fields are classified by a choice of isotropy genera- 
tor and SO{5) is distinguished as the simplest model with a non-Abelian resid- 

^■f^ , ual (little) group, SU (2) X U (1), and which admits globally regular particle-like 

■ solutions. We exhibit an algebraic gauge condition which normalises the resid- 
ual gauge freedom to a finite number of discrete symmetries. This generalises 
the well-known reduction to the real magnetic potential w{r,t) in the original 

CJ '. SU (2) YM model. Reformulating using gauge invariant polynomials dramat- 

ically simplifies the system and makes numerical search techniques feasible. 
> * ^ ' We find three families of embedded SU{2) EYM equations within the 50(5) 

{2J[j| system, one of which was first detected only within the gauge-invariant polyno- 

mial reduced system. Numerical solutions representing mixtures of the three 
SU{2) sub-systems are found, classified by a pair of positive integers. 

I ■ 
> 

in 

G^ ■ 1. Introduction 

m 

I I The LHC experiments currently underway at CERN are expected to settle promi- 

■ nent vexing questions such as the origin of the Higgs mechanism, the existence of 
^\ I supersymmetry in nature, and the meaning of dark matter. However, regardless of 

i the discoveries still to flow from the LHC, we can with confidence predict that the 

governing equations will incorporate the Einstein equations for the gravitational 
field and the Yang-Mills equations for some non-Abelian gauge group. From this 
perspective, the spherically symmetric solutions of the S'[/(2)-EYM equations dis- 
covered numerically in [110 are particularly important, for several reasons. First, 
the particle- like properties exhibited by the solutions, namely static, asymptotically 
flat, and globally regular with spatial topology M.^ , confound previous expectations, 
based on the known non-existence results for the vacuum Einstein equations and 
for the YM equations separately, that such solutions could not exist. On the other 
hand, they confirm Wheeler's "geon" hypothesis [15], that localised semi-bound 
gravitational (and with hindsight, Yang-Mills) configurations might exist. 

Numerical and perturbation results [17,18] show that the SU{2) EYM spherically 
symmetric solutions may be viewed as an unstable balance between a dispersive YM 
and the attractive gravitational force. However, there are many aspects of these 
systems which have not yet been studied in depth, and it seems premature to 
conclude that instability is inevitable — see the review [16]. 
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"'^ Existence of global solutions was proved rigorously in [12, 13]. 
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The most detailed results have generally been obtained in the static spherically 
symmetric setting with the simplest YM gauge group SU{2), which exhibits many 
of the basic properties of the general non-Abelian group models. The spherically 
symmetric reduction of the EYM equations can be viewed as a 2D EYM-Higgs 
system with a residual gauge group, G^^ , a residual Higgs field (to be defined below), 
and a "Mexican hat" potential. Most attention has focused on the case where the 
gauge group is SU (n) and the residual group is Abelian, primarily because in this 
case it is possible to completely fix the gauge freedom. If the residual group is 
non-Abelian, then it is known [4, 5] that the issue of gauge fixing becomes much 
more challenging. 

To better understand gauge fixing and solutions to the the static spherically 
symmetric EYM equations we study G = SO{5), which is the simplest gauge group 
that supports a globally regular model with non-Abelian residual group. This model 
was discussed in [4] , but the questions of gauge fixing and the existence of solutions 
were not resolved. We solve this problem here by using a related system satisfied by 
polynomials in the gauge fields which are invariant under the action of the residual 
gauge group. 

2. Static spherically symmetric field equations 

Spherical symmetry for Yang-Mills fields is complicated to define because there 
are many ways to lift an iS'0(3) action on space-time to an action on the Yang- 
Mills connections. For real compact semisimple gauge groups G, it was shown 
in [1,3] that equivalent spherically symmetric Yang-Mills connections correspond 
to conjugacy classes of homomorphisms of the isotropy subgroup, U{1), into G. The 
latter, in turn, are given by their generator A3, the image of the basis vector T3 of 
su(2) (where {r^, z = 1, . . . , 3} is a standard basis with [n, r,] = e'^^jTk), lying in an 
integral lattice / of a Cartan subalgebra f)o of the Lie algebra go of G. This vector 
A3, when nontrivial, then characterizes up to conjugacy an su(2) subalgebra0 

With the vector A3 e go fixed, a gauge and a coordinate system (t, r, 6, (f)) can 
be chosen [1,3] so that the metric and gauge potential take the form 

ds^ = -S{r f - dt^ ' dr^ + r^ {d0^ + sin^ d<i>^) , 

ancQ 

A = Ki [r)de + (A2 (r) sin 6* -f A3 cos 9)d(j), 
where Ai(r), A2(r) arc go"Valucd maps satisfying 

(2.1) [A3,Ai(r)] A2(r), and [A2(r), A3] - Ai(r). 
It is more convenient to use the following variables, 

(2.2) Ao := 2iA3, and A± := ^pAi - iA2, 

which lie in g, the complcxification of go. The conditions (|2.ip on Ai, A2 then 
imply that the connection functions h.± are valued in the residual Higgs bundles, 
V±2, defined by 

(2.3) A±(r) e y±2 := { ^ € 9 I [Ao, X] = ±2X} . 

^This ignores some interesting effects due to the fact that 50(3) is not simply connected. For 
an analysis of 50(3) actions on SU{2) bundles see [2]. 

■^Here we are assuming the ansatz, that the gauge potential is "purely magnetic". 
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The residual (gauge) group G^^ is defined as the connected Lie subgroup with 
Lie algebra 

/■^ :={Xe0o|[A3,X] = O}. 
Wc call the residual group (non-)Abelian if this is a (non-)Abelian Lie algebra. For 
g^^ to be non-Abelian it is necessary but not sufficient that Aq lies on the boundary 
of a Weyl chamber of f), a Cartan subalgebra of g. 

In the variables (|2.2p . the EYM equations become [6] 

(2.4) m' = ^||A;|p + ^||Ao-[A+,A_]|p, 

(2.5) r2(5W+)' + ^(^A+-i[[A+,A_],A+]) =0, 

(2.6) [AV,A_] + [A'_,A+] =0, 

(2.7) S-^S' = -\\A'\\^, 

r 

where iV = (1 — 2m/r). The S parameter can be decoupled from the system, 
decreasing the order of the overall system by one but introducing a first order term 
in (|2.5p . Also note that the second term in l|2.5p is proportional to the gradient of the 
second term in (|2.4p . Requiring that the solutions are regular and asymptotically 
flat gives boundary conditions 

[A+,A_]=Ao, and [[A+, A_], A+] = 2A+ 

at r = 0, and as r ^ oo, respectively. 

The 1 1 • I |-norm is proportional, on each irreducible component of g, to the the real 
part of the Hermitean inner-product derived from the Killing form [6] . Multiplying 
II • II by a constant factor leads via a global rescaling of m and r to the original 
equations. Specifically, we have: 

Proposition 1. If {m{r),h+{r)) satisfies (|2.4p - (|2.6p . then {am{r / a), h+{r / a)) 
satisfies the equations obtained by replacing || • || in equations (|2.4p - (|2.6p with a|| • ||. 

Proof. Substitution. □ 

The EYM system (|2.4p - (|2.6p typically admits subfamilies of solutions which also 
satisfy the original SU{2) but with rescaled norm || • || (see Section 5). We may 
regard such solutions are equivalent to the original SU(2), since all values are 
consistently rescaled from the values in the original reports [7-9] . 

3. An 50(5) model 

We now specialize to the gauge group G = SO {5). The complexified Lie algebra 
g = so(5, C) has the Cartan decomposition 

g = (je0Ce„, 

where f) = spauj^l/iaj , /iq^] is the Cartan subalgebra, R — {a±i,i = 1, ... ,4} is a 
root system in f)* with root diagram as in Figure [1] 

The reduced spherically symmetric EYM model depends strongly on the choice 
of isotropy generator Aq. For G = SO{5), we choose [4] 

Ao = 2Hi. 
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Figure 1. Root Diagram for so(5,C) 



With this choice, the residual group is SU{2) x U{1), while (see 
and hence 

(3.1) A+(f) = wi{r)eai + W2{r)ea2 + W3ir)eas, 

(3.2) ^-(r) =Wi(r)6o_i +W2ir)ea_2 +W3{r)ea_3. 

Substituting (|3.ip and (|3.2p into (|2.4p - (|2.6p . we find that the static spherically 
symmetric EYM equations in terms of the Wi are 

(3.3) ,ri' = ^i\w[\' + \w'2\' + \w',\') + ^V, 

(3.4) r^Nw'l + {2m - -V)w[ + wi{l - \wi\^) + W2{^^ 

r 2 

(3.5) r^Nw'^ + (2m - -r)w'^ + ^2(1 - - kiP - ksH + 

r 2 



- W1W2) 


= 0, 




= 0, 


- W3?U2) 


= 0, 




-0, 


- W3UJ3) 


= 0, 


- W2W3) 


= 0, 



(3.6) r^Nw'-i + {2m^-V)w'3 + W3{l-\w3\^)+W2i- ^ 

(3.7) w'lWi — WiW']^ + ?«2^2 ~ W2U'2 + W3' 

(3.8) w[wi — wiw'i — {w'^u 

(3.9) w'iW2 — 'wiw'2 — {w'2U 

where 

V ii^^ + ii^^ + |.2P(krP + hap + ^ - 1) - 

= i + 1^2^ + kaP - 2) V i {\wi\^ - Iwap)^ + 1 \w1W2 - W2W3\'^ 
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The boundary conditions are 

(3.10) |u.i|2 + |u;2|' + h3p = 2, 

(3.11) \wi\-\w3\=0, 

(3.12) W1W2 — W2W3 = 0, 



at r = 0, and 

(3.13) wiil - \wi\'') + W2{^^ - = 0, 

(3.14) W2il - ^^) - W2{\wi\^ + \W3\'')+W2WIW3 = 0, 

(3.15) W3{1 ~ IwaP) + W2 C"'^^^ ~ W3W2) ^ 0, 
as r ^ 00. 

The sets of boundary conditions (|3.10p - (|3.12p and (|3.13[) - p.l5[) define algebraic 
varieties on and as such they do not restrict the initial and final values to one 
particular point in (or indeed a discrete set of such points) . Given a point lying 
on one of these varieties there will in general be a set of gauge-equivalent points but 
it is not immediately obvious what these equivalence classes of points are from the 
defining equations. This is a problematic feature of all models with a non-Abelian 
residual gauge group and will be discussed elsewhere. 

4. Gauge fixing and the reduced equations 

For the SO{5) model, a direct computation shows that the ansatz 

(4.1) A+(r) = w(r)eQi + 06^2 + w(r)ea3, u(r), i;(r) e M 

automatically satisfies the constraint equations ()2.6p (equivalently (|3.7p - (|3.9p ). while 
the equations (iMl-jMl) (equivalently (ISSll-lISH)) reduce to 



(4.2) 



r^Nu" + (^2m - (1^^0_Lii^^^ ^' + ^(1 _ ^ 

(4.3) r^Nv" + (^2m - ~ + ~ ^ v' + v{l - v^) = 0, 

(4.4) m' = ^{{ur + {vT) + ^((1 - + (1 - ^'')')' 
with boundary conditions 

(4.5) u^+v'^ = 2, \u\ = \v\ 
at r = and 

(4.6) m(1-m2) = o, v{1-v^)=0 
as r — > oo. 

The above results show that (|4.ip is a consistent ansatz for the static spherically 
symmetric EYM equations. In fact the equations that result could have been ob- 
tained from a model with gauge group SU{2) x SU{2). However, it is not obvious 
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that the choice (|4.1|) is equivalent to fixing a section of the gauge fields. To see this 
we use the fact that the following two polynomials 

K, :==2||A+||2=2(|ii;i|2 + |t.2|2 + |y;3|2), 

K2 2||[A_, A+]||2 =4|ii;i|4 + 2\w2\^ + A\w:i\^ 

+S\wi\^\w2\^ + 8|w3nw2p - 8Re[w;iwW], 
are a complete set of generator^ for the ring of residual-group- invariant polynomials 
in A_|_. When {wi, W2, W3) = (m, 0, v) for real it, v, these become 

Ki=:2{u^ + v'^), and K2 = '0^' + v'^). 
These equations can be inverted to give 



u = ±,\^—±,y j , and v^±2\^—T3^ ^ J 

which are well-defined for all Ki , K2 since the inequalities 

2K2 - Kf > 0, and K2 - Kl < 

follow directly from the definition of Ki and K2 ■ The ±-signs amount to choosing 
which function is u and which is v and deciding on an arbitrary convention for their 
initial values. 

5. Numerical results 



From the form of the reduced equations (|4.2[) - (|4.4[) . it is clear that we can make 
two distinct simplifying assumptions u ^ v or v = ±1 (equivalently u = ±1)- 

(i) Setting u = v, the reduced equations (|4.2p - (|4.4p become 

r^Nu" + (2m ~ ii— u' + u{l - u^) = 0, 

2r^ 

which is the SU (2) equation and has the well-known Bartnik-McKinnon family of 
solutions [7,11]. We note that these correspond to the embedded SU{2) solutions as 
they satisfy A±(r) = u(r)Q.± where il± are fixed vectors satisfying [Aq, ft±] = 2n± 
and [ri+, ri_] = Aq. 

(ii) Setting v ~ 1, the reduced equations (I4.2p - (|4.4p become 

r^Nu" + (^2m - + "(1 ' - 0, 

As in Proposition [1] rescaling these equations by a/2 in m and r results in the 
S'?7(2)-equations above, giving another family of solutions where the u field is a 
radially-scaled Bartnik-McKinnon solution with smaller mass and the v field re- 
mains constant. 



"This can be established by computing the Molicn function to determine the appropriate orders 
of the polynomials and then verifying Ki and K2 are algebraically independent. See [10] for some 
useful formulae. 
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12 3 



Figure 2. A plot of the solutions in the space of the invariant 
polynomials. All solutions start at (4, 8). The type (i) SU{2) so- 
lutions lie on the K2 = ^Kf curve, while the type (ii) SU{2) 
solutions he on the K2 = {Ki — 2)^+4 curve. The type (iii) so- 
lutions are curves in the combined interior of the three parabolas 
touching the boundaries tangentially (an illustration of the (2,1) 
type (iii) solution is shown) 



For the above two special situations, existence of a countably infinite number 
of solutions is guaranteed by the existence theorems in [7, 12, 13]. Integrating the 
reduced equations (|4.2p - (|4.4p numerically, we also find evidence for solutions that 
are not of the type (i) or (ii) above. We will call these solutions type (iii). Figure [2] 
indicates where these solutions lie in the space of invariant polynomials. All three 
types of solutions can be characterized by the number of nodes of the functions u 
and V. Figure [3] enumerates all of the solutions that were found numerically or that 
are known to exist analytically. The type (i) solutions lie on the diagonal n„ = riy, 
while the type (ii) solutions lie on either the horizontal n„ or vertical axis 
riu = 0. The remaining solutions are of type (iii) and are indicated on the node 
diagram by the larger circles. We conjecture, based on the diagram, that for each 
(n„,n^) G Z>o X Z>o there exists a solution to the reduced equations (|4.2[) - (|4.4p 
satisfying the boundary conditions (|4.5p - (|4.6p . 

For the numerical analysis we use the "shooting to a fitting point" method [14], 
using the Matlab routine ode45 which is a fourth-order Runge-Kutta solver. Using 
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1 2 3 4 5 

Figure 3. Node Diagram for solutions to the reduced equation 



the power series expansions 

9 9 3 

u = l + air^ + (^aj + ^aiaj + —ajy^ + 



^5 
2 

1 + a2r2 + Ga2al 
5 



10 
3 

10 



4y 



0(r6), 



{al+al)r^ + Qal + ^4y + 0{r'), 
5 5 



near r = 0, and 



±1 



h 

r 

&2 



0(4), 



« = ±1 + 7 + 0(^:2)' 

m = TOoo + O(^), 

near r = 00. We shoot from the origin out trying to maximise (by searching in the 
parameters ai evaluated at r = 0.01) the r-value that is attained before the solution 
violates one of the necessary conditions for a global smooth solution. Suitable 
large r-values will then give an indication what the hi and moo parameters should 
approximately be for the shooting back from some large r-value (we use r = 10000). 
We then use these parameter values as an initial guess in the "shooting to a fitting 
point" search, where the fitting was done at r = 10. 



6. Discussion 

We have found three solutions using the numerical technique from the previous 
section. When combined with the two families of SU{2) type solutions, they suggest 
that there is a solution for each point {Tiu,njj) on the nonnegative integer lattice, 
where n„ , Uy denote the number of nodes of u and v respectively. The numerical 
values arc given in Table [1] and some type (iii) solutions are displayed in Figures [4] 
andEl 
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It is only with the reduced variables we have presented that these solutions could 
be found. As an indication of the difficulties otherwise, consider the particulars of 
the shooting method technique. In the original lu-variables the ODE system is 
given in terms of four complex variables and the mass, m. This leads to eight 
real shooting parameters near r = and nine shooting parameters near r = oo. 
Furthermore, after finding a promising approximate solution by shooting outward 
from r = and another by shooting inward from r = oo, they will generally need to 
be gauge rotated to match at the fitting point (if indeed they are gauge equivalent- 
solutions). This all amounts to a difficult optimization problem in 17 variables. 
This problem was discussed in [4] where some approximate solutions were obtained 
by shooting but a global solution was not found. 

In contrast, shooting for the reduced variables has only five parameters to search 
in, two parameters near r = and three near r = oo with no need to gauge 
rotate possible solutions at the fitting point. The removal of degeneracy from the 
numerical problem and minimization of the number of shooting parameters is what 
has made it tractable to find the solutions presented here. 
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-0.65173 
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-8.8639 
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-1 


1 


-1.246544 


-0.155699 


42.52411 


-13.88445 


3 
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-1 


-1 


-0.69704 


-0.69704 


58.9326 


58.9326 
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(i) 


0.99924 


1 


1 


-0.70488 


-0.70488 


-366.335 


-366.335 


4 
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0.99988 


-1 


-1 


-0.70617 


-0.70617 
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2251.89 
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0.99998 
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-0.70638 


-0.70638 


-13817.4 


-13817.4 
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Table 1 . Table of Numerical Parameters 



-2-101234 

Figure 4. The (2, 1) type {in) solution. The mass, to, increases 
monotonicallv from to 0.924. u starts and ends at 1 with two 
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Figure 5. The (3,2) type (Hi) solution. The mass, to, increases 
monotonically from to 0.987, u goes from 1 to —1 with three 
nodes and v goes from +1 to +1 with two nodes. 
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